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Quantum optics has been a major driving force behind the rapid experimental developments 
that have led from the first laser cooling schemes to the Bose- Einstein condensation (BEC) of 
dilute atomic and molecular gases. Not only has it provided experimentalists with the necessary 
tools to create ultra-cold atomic systems, but it has also provided theorists with a formalism and 
framework to describe them: many effects now being studied in quantum-degenerate atomic and 
molecular systems find a very natural explanation in a quantum optics picture. This article briefly 
reviews three such examples that find their direct inspiration in the trailblazing work carried out 
over the years by Herbert Walther, one of the true giants of that field. Specifically, we use an 
analogy with the micromaser to analyze ultra-cold molecules in a double-well potential; study 
the formation and dissociation dynamics of molecules using the passage time statistics familiar 
from superradiance and superfluorescence studies; and show how molecules can be used to probe 
higher-order correlations in ultra-cold atomic gases, in particular bunching and antibunching. 



I. INTRODUCTION 



Quantum optics plays a central role in the physics of quantum-degenerate atoms and molecules. Laser light 
and its coherent and incoherent interactions with atoms are ubiquitous in t hese experiments, and the tools that 



have culminated i n the achievement of Bose-Einstein condensation (BEC) l|Anderson et al 



1995; 



Davis et al 



1995; 



Bradley et al 



1995|) were first studied and understood in quantum optics. Indeed, the deep connection between 
quantum optics and cold atom physics was realized well before the first experimental realizations of BEC, both at the 
experimental and theoretical levels. On the theory side, there are (at least) two important reasons why quantum optics 
methods are well suited for the study of cold atoms systems. First, bosonic fields have direct analogs in electromagnetic 
fields, which have been extensively studied in quantum optics. Second, for fermions the Pauli Exclusion Principle 
restricts the occupation of a given mode to zero or one, and these two states — mode occupied or empty — can often 
be mapped onto a two-level system, as we shall see. As a result, ma ny situations familiar from quantum optics are 



also found in col d-atom systems, including matter-wave interference And^w^^t^al 



wave amplifiers Ijlnouve et al 



1999b; 



Ketterle and H.-J. Miesner 



matter-wave beam splitters IjBurgbacher and Audretsc. 



ch, 



1997 



Kozuma et al 



1999), four- wave mixing IjChrist et al 



9971) atom lasers and matter 



1999 



Law and Bigelow 



200.i 



1998) 



Lenz et al 



2 



1993; 



Meiser et al 



2005a: 



Dicke superradiance l|Inouve et al 



1999e 



2003; 



Moore et a 



Moore and Mevstre , 



1999; 

3, 



Roio et al 



1999; 



Search et al 



2002b), and 



19901) . to name a few. 

At the same time the physics of ultra-co l d atoms is much richer than its quan tum-optical counterpart since atoms 



can be either fermions (DcMarco and Jin 



1998, 



1999; 



Hadzibabic et al 



2003) or bosons and have a rich internal 



structure. In addition, the in teraction between atoms can be tuned relatively easily on fast time scales using; for 



instance Feshbach resonances (Duine and Stool 



Timmermans et al 



2004; 



Purr et al 



2004; 



1999|) or two-photon Raman transitions IjTheis et al 



noirveet^al 



2004; 



1998, 



2004; 



Wvnar et al 



Stan et al 



2004: 



2000). Indeed, some 



of the most exciting recent developments in the physics of ultra- cold atoms are related to the coher ent coupling 



of atoms to ultra-cold molecules by means of Feshba ch resonances ( Diir r et al 



association (Kcrman ct al 



2004: 



Wvnar et al 



Diirr et al. 


2004 


Rceal et al. 


2003) 



2000|) . Both bo sons and fermions have been successfully converted into 



molecules. In both cases BEC of molecules has been observed l|Donlev et al 



2003; 



Zwierlein et al 



2002; 



Greiner et al 



2003; 



Jochim et al 



2003) , and the long-st anding question of the BEC-BCS crossover is 



mentally and theoretically in those systems 



Regal et al 



2004a: 



Timmermans et al 



2001 



I Bartenstein ct a 



Zwierlein et al 



2004; 



Holland et al 



2001 



being investigated 



experi- 



Ohashi and G riffin , 



2002; 



2004). Other developments with close connections with 



2002a b; 



2004; 



Jaksch 



Walthe 



et al 



1998|) . which 



1992) and leads in 



quantum optics include the trapping of atoms in optical lattices {Gremgr^^d^ 
play a role closely related to a high-Q resonator in cavity QED ( Se arch et al. 
addition to fascinating connections with condensed matter physics and quantum information science. 

With so many close connections between the physics of quantum-degenerate atomic and molecular systems and 
quantum optics, it is natural and wise to go back to the masters of that field to find inspiration and guidance, 
and this is why Herbert Walthcr's intellectual imprint remains so important. This brief review illustrates this point 
with three examples. Section ^ shows that the conversion of pairs of fermions into molecules in a double-well 
potential can be described by a generalized Jaynes-Cummings model. Using this equivalence, we show that the 
dynamics of the molecular field at each site can be mapped to that of a micromaser, one of Herbert Walther's most 



remarkable contributions (Mcschcde et al 



1985). Section IlIII further expands on the mapping of ultra-cold fermion 
pairs onto two-level atoms to study the role of fluctuations in the association and dissociation rates of ultra-cold 
molecules. We show that this system is closely related to Dicke superradiance, and with this analogy as a guide, we 
discuss how the passage time fluctuations depend sensitively on the initial state of the system. In a third example 



inspired by Herbert Walther's work on photon statistics and antibunching IfBrattke et al 



Rempe et al 



1990; 



Rempe and Walther 



2001; 



Krause et al 



1989 



1990) scction lTVl analvzes how the statistics of their constituent atoms affects 



3 

the counting statistics of molecules formed by photo-association. We compare the three cases where the molecules 
are formed from a BEC, an ultra-cold Fermi gas and a Fermi system with a superfiuid component. The concept of 
quantum coherence developed by R. J. Glauber and exploited in many situations by H. Walther and his coworkers, 
in particular in their studies of resonance fluorescence, are now applied to characterizing the statistical properties of 
the coupled atom- molecule system. Finally, section further elaborates on these ideas to probe spatial correlations 
and coherent properties of atomic samples, and we find that the momentum distribution of the molecules contains 
detailed information about the second-order correlations of the initial atomic gas. 



II. MOLECULAR MICROMASER 



Ultra-cold atoms and molecules trapped in optical lattices provide an exciting new tool to study a variety of 
physics problems. In particular, they provide remarkable connections with the condensed matter of strongly cor- 
related systems and with quantum i nformation science, a very well controlled environment to study processes such 



as photo-association l|Rvu et al 



2005|) . and, from a point-of-view more directly related to quantum optics, can be 
thought of as matter-wave analog of photons trapped in high-Q cavities. In particular, the high degree of real-time 
control of the system parameters offers the opportunity to directly experimentally study some of the long-standing 
questions of cond ensed matter phys i cs, su ch as the ground state structure of certain models and many-body dy- 



namic properties iJaksch and Z ollcr 



2005). The coherent formation of molecules in an optical lat tice via either 



Feshbach resonances and two-photon Raman photo-association has been studied both theor etically ( Da mski et al 



2003; 



tally 



Esslingerj 



Kohl et al 



i Rvu et al 



d Mo 



200.4 



mer 



2003; 



Rom et al 



Jaksch et al 



2004; 



Rvu et al 



2002: iMo 
al, 200.4 



mer 



2003; 



Moore and Sadeehpour 



Stoferle et al 



)5) 



2003) and experimen- 



2005). In particular, the experiment of Ref. 



2005) observed reversible and coherent Rabi oscilla tions in a gas of coupled atoms and molecules 



The idea of the molecular micromaser (Search et al 



2003) relies on the observation that, as a consequence of 



Fermi statistics, the photo-association of fermionic atoms into bosonic molecules can be mapped onto a generalized 
Jaynes-Cummings model. This analogy allows one to immediately translate many of the results that have been 
obtained for the Jaynes-Cummings model to atom-molecule systems. In addition, the molecular system possesses 
several properties that have no co unterpart in the quantum optics analog , giving rise to interesting generalizations of 



the original micromaser problem (Filipo wicz et al 



1986; 



Guzman et al 



1989 



Meschede et al 



1985 



Remoe et al 



1990). One of these new features is the inter-site tunneling of atoms and molec ules between adjace nt lattice sites, 



leading to a system that can be thought of as an array of molecular micromasers IjSearch et al 



2004). 



To see how this works, rather than treating a full lattice potential we consider the dynamics of the molecular field 
in the simpler model of a coupled atom-molecule system in a double-well potential. We first show that inter-well 
tunneling enhances number fluctuations and eliminates trapping states in a manner similar to thermal fluctuations. 
We also examine the buildup of the relative phase between the two molecular states localized at the two wells due 
to the combined effect of inter-well tunneling and two-body collisions. We identify three regimes, characterized by 
different orders of magnitude of the ratio of the two-body collision strength to the inter-well tunneling coupling. The 
crossover of the non-equilibrium steady state from a phase-coherent regime t o a phase-incoh erent regime is closely 



related to the phase locking of condensates in Josephson-type configurations ([Leggett , 



2001), while we consider an 



open quantum system with incoherent pump and molecular loss which results in a dissipative steady state. 
A. Model 

We consider a mixture of two hyperfine spin states \a =|, |) of fermionic atoms of mass nif trapped in a double- well 
potential at temperature T = 0, which can be coherently combined into bosonic molecules of mass nib via two-photon 
Raman photo-association. If the band-gap of the lattice potential is much larger than any other energy scale in the 
system, the fermions and molecules occupy only the lowest energy level of each well and the number of fcrmions of a 
given spin state is at most one in each well. 

In the tight binding approximation, the effective Hamiltonian describing the coupled atom-molecule system is 

i— l,r 

where 

H oi = h(uj b + 5)hi + frujf(h^ l + n^) + ^hUbh^fii - 1) 

+ hU x hi(n^i + nii) + hUfh^nii, (2) 
Hit = hx(t)b\c^c li + H.c, (3) 
Ht = —hJbb\b r — hJf(cyC<i r + cj ; C| r ) + H.c. (4) 

Here c a i and hi, i = l,r, are the annihilation operators of fermionic atoms and bosonic molecules in the left (I) and 
right (r) wells, respectively. The corresponding number operators hi — b\bi and n a i = c^c CT j have eigenvalues rii and 
n a i, respectively, and hu>b and hujf are the energies of the molecules and atoms in the isolated wells. 

The terms proportional to Ub, U x , and Uf in H$i describe on-site two-body interactions between molecules, be- 
tween atoms and molecules, and between atoms, respectively. The interaction Hamiltonian Hn describes the con- 



version of atoms into molecules via two-photon Raman photo-association. The photo-association coupling constant 
x(t) is proportional to the far off-resonant two-photon Rabi frequency associated with two nearly co-propagating 



lasers (Heinzen et al 



2000), and S is the two-photon detuning between the lasers and the energy difference of the 
atom pairs and the molecules. The tunneling between two wells is described by the parameters Jb and Jf in the 
tunneling Hamiltonian Ht- 

The molecular field is "pumped" by a train of short photo-association pulses of duration r, separated by long 
intervals r»r during which the molecules are subject only to two-body collisions and quantum tunneling between 
the potential wells, as well as to losses due mainly to three-body collisions and collisional relaxation to low- lying 
vibrational states. In the absence of inter-well tunneling, this separation of time scales leads to a situation very 
similar to that encountered in the description of traditional micromasers, with the transit of individual two-level 
atoms through the micromaser cavity replaced by the train of photo-association pulses. 

The dynamics of the molecular field in the double-well system is governed by the following four mechanisms: 
(i) Coherent pumping by injection of pairs of fermionic atoms inside the double-well. This process is the analog of 
the injection of two-level atoms into a micromaser cavity. The injection of pairs of fermionic atoms into the doub l e- well 



potential can be acc omplished e.g. by Raman transfer of atoms from an untrapped internal state ( Jaksch^eta 



Mandel et al 



1998; 



2002a). 



2004). This results in the pumping of fermions into the double well at a rate T l|Search et al 
We assume that for times T ^> r _1 , a pair of fermions has been transferred to the two wells with unit probability. 



that is, the state of the trapped fermions in well i is 



(5) 



(ii) Molecular damping, which is the analog of cavity dampi ng. During the ti me intervals T when the photo 



2003). The decay of the molecules is 



association lasers are off, the molecular field decays at rate 7 ({Search et al 
due to Rayleigh scattering from the intermediate molecular excited state, three-body inelastic collisions between a 
molecule and two fermions, and collisional relaxation from a vibrationally e xcited molecular state to deeply bound 



states. These loss mechanisms c an be modeled by a master equation, see e.g. IjMevstre and III , 



2004: 



Scully and Zubairv 



1997) 



1999; 



Mivakawa et al 



(iii) The application of a train of photo-association pulses. This mechanism is formally analogous to the Jaynes- 
Cummings interaction between the single-mode field and a sequence of two-level atoms traveling through the microwave 
cavity in the conventional micromaser. As already mentioned, we assume that these are square pulses of duration 



r and period T + r, with r much shorter than all other time scales in this model, r <C J b *,7 1 . This assumption 



is essential if we are to neglect damping and tunneling while the photo-association fields are on. The change in the 
molecular field resulting from atom-molecule conversion is given by 

Fi(T)p b = Tr a [Ui{T)Mt)uHr)], (6) 

where p ab is the total density operator of the atom-molecule system and Tr a [ ] denotes the trace over the atomic 
variables, Ui(r) = exp (— i/ijT / K) being the evolution operator for a single-well Hamiltonian, 

hi = Hgi + Hh- 

The key observation that allows us here and below to bui ld a bridge from the cold atoms and molecular system to 



quantum optics systems is that by means of the mapping (Anderson, 



1958) 



°+i = c \Ai> 

the atomic degrees of freedom take the form of a fictitious two level system. The operators cr + i,<7_i, and a Z i can 
be interpreted as the raising and lowering operators of the fictitious two-level atom and the population difference, 
respectively. The single-well Hamiltonian takes the form 

hi = h(uj b + U x ) hi + h(uo f + U x hi) a zl 

+ h (x{t)b\a-i + X *(t)kv+i) + ^U b ni(hi - 1) (8) 

where we have dropped constant terms and we have redefined u>b and ujf according to u>b + 5 — > LOb and ujf + U f /2 — > ujf. 

The Hamiltonian hi is Jaynes-Cummings-like, and for \ = const., the resulting dynamics can be determined within 
the two-state manifolds of each well {|ei,nj), \gi,rii + 1)} by a simple extension of the familiar solution to the Jaynes- 
Cummings model. Within each manifold the system undergoes Rabi oscillations. 

Since tunneling is neglected during the photo-association steps, the two wells are independent of each other and 
identical to each other. The resultant molecular gain is then modeled by independent coarse-grained equations of 
motion for the reduced density matrices of each molecular mode. 

(iv) The unitary time evolution of the molecular field under the influence of two-body collisions and quantum 
tunneling, a process absent in conventional micromasers. During the intervals T it is governed by 

- w = — n [H b ,pb], (9) 
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where the Hamiltonian 

H b = -hJ b {b\l r + blk) + h^j-(ni - h r f (10) 

contains tunneling and collisions. In Eq. I|10|) . we have neglected terms that are functions only of N = hi + n r , a step 
justified as long as the initial density matrix is diagonal in the total number of molecules in the two wells. 

Combining the coherent and incoherent processes (i) to (iv), we obtain the full evolution of the molecular field 

Ufo + \ £ Mr) h] p b -\ [H b , p b ] , (11) 

l,r l 7 r 

where p b is the reduced density matrix of the molecules. The initial condition for the molecules is taken to be the 
vacuum state. Because the molecular pumping and decay is the same in both wells, the density matrix p remains 
diagonal in the total number of molecules in the two wells for all times. This is a generalization of the micromaser result 
that the photon density matrix will remain diagonal if it is initially diagonal in a number state basis jFilipowicz et al. , 
1983). 



B. Results 

The master equation describing the molecular micromaser dynamics contains six independent parameters: the 
number of photo-association cycles per lifetime of the molecule, N ex = I/7T; the "pump parameter" = \J N ex \y\r\ 
the two-body collision strength and tunneling coupling strength per decay rate, u b = U b /^i and tj — J b /"f] and finally, 
the detuning parameter 77 = (2ujf — uj b )/2\x\ and the nonlinear detuning parameter (3 = (2U X — J7(,)/2|x|- 

In our model, the atomic and molecular level separations in the wells are required to be much larger than the 
relevant interaction energies, 



fkj b > U b (ni)((ni} - 1), \x\VJfk), 



(12) 



hj) being the 



I Greiner et al 



average number o f molecules in well i. A c omparison with actual experimental parameters 



2002a: 



Jaksch et al 



1998: 



Mivakawa et al 



2004) shows that these conditions are satisfied as long 



as the number of molecules doesn't exceed 10. In addition, the neglect of inter-well tunneling and damping effects 
during the photo-association pulses requires that 



(13) 



This condition is satisfied in typical experiments. 
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In the remainder of this section, we discuss the dynamics of the molecular field obtained by direct numerical 
integration of the master equation with a Runge-Kutta algorithm until a dissipative steady state is reached. For 
simplicity, we confine our discussion to exact resonance only, 77 = j3 = 0, and a fixed value of N ex = 10. 

1. single-well molecular statistics 

We first discuss the statistics of a single- well molecular mode, which is given by tracing over the full density matrix 
with respect to degrees of freedom of the other localized mode as 



We note that off-diagonal elements of the density matrix for a single well are zero. Since the initial state of the 
molecules in each well is the same, i.e. the vacuum state, and Hb is invariant with respect to the interchange I «-> r, 
the molecule statistics for left and right wells are identical, P(ni) = P{n r ) = P(n). 

Figure^shows the steady-state average number (hi) is plotted as a function of the pump parameter for ui, = 0, 
N ex = 10. In the absence of inter- well tunneling, corresponding to Fig.^a), the result reproduces that of conventional 
micromasers, with a "lasing" threshold behavior at around f=a 1 and an abrupt jump to a higher mean occupation 
at about the first transition point, ~ 2n. The former effect is not affected by the tunneling coupling as shown in 
Fig. ^b). However, the latter abrupt jump disappears in the presence of interwcll tunneling. This is because the 
coupling to the other well leads to fluctuations in the number of molecules in each well and has an effect similar to 
thermal fluctuations in the traditional micromaser theory. The enhancement of fluctuations can also be seen in Fig. [21 
where the Mandel (^-parameter 




is plotted as a function of O. 

It is known that in the usual micromaser the sharp resonance- like dips in (hi) and Q are attributable to trapping 
states, which are characterized by a sharp photon number. For the specific value of = y/Sn, as shown in Fig, ffla). 
the number probability does not reach beyond number state |rij = 1) in the case of tj = 0. As shown in Fig. Ufb), 
the tunneling coupling makes possible transitions into higher number states and eliminates the trapping state in a 
manner similar to thermal fluctuations in the conventional micro maser. 




(14) 
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d 
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FIG. 1 (fii) versus Q/n for Ub = and iV e:! ; = 10, and for (a) t j = and (b) £j = 5. 



2. phase coherence between two micromasers with tunneling coupling 



So far we have discussed the single-well molecule statistics and how it is affected by inter-well tunneling. Now 
we turn to a more detailed discussion of the phase coherence between the two localized modes. It is very useful to 
divide the par ameter space o f the ratio of the two-body collision strength to the inter-well tunneling coupling into 



three regimes IjLeggett . 



20011) : "Rabi-regime" (u b /tj < (TV) -1 ); "Josephson-regime" {(N)- 1 < u b /tj < (N)); and 
"Fock-regime" ((N) <C Uh/tj); where (N) denotes the average total molecule number. 
The analysis of the relative coherence of the molecular fields in the two wells is most conveniently discussed in terms 



of the angular momentum representation 



J+ 
J_ 

Jz 

J 2 



J x - ij v = b\bi, 



1 



(blk - btbr), 



2 

N ( N 



(15) 
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FIG. 4 {Jx)/{nj) versus Ut/tj for O = n and tj = 2.5. 

The symmetry of the density matrix with respect to the two wells furthermore implies that (J z ) = ( J y ) — 0. The 
first-order coherence between the molecular fields in the left and right potential wells is then given by (J x ). Figured 
shows the normalized steady-state first-order coherence (Jx}/(nj) as a function of Ub/tj for 9 = it and tj = 2.5. 
(J x ) is suppressed in both the Rabi and Fock regimes and has an extremum at \ub\/tj ~ 0.6. In the Fock regime, 
\ub\/tj 3> (N) ~ 10, the nonlinearity in Hb dominates and reduces the coherence between the localized states of 
each well. We note that the average occupation numbers for each well are relatively unaffected by Ub/tj, with 
(hj) = (N) /2 = 4.78 - 4.87 for \u b \/tj = 10 2 - 10~ 2 - 5 . 

The reason why the first-order coherence is suppressed in the weak coupling limit, Ub — 0, can be understood 
as follows. The expectation value (J x ) corresponds to the difference in occupation numbers between the in-phase, 
b s = (pi + b r )/y/2, and out-of-phase, b a = (bi — 6 r )/\/2, states of the localized states of each well, J x = b\b s — b\b a . 
Since the bandwidth of the photo-association pulse is larger than their energy splitting, 1/t» those states are 
equally populated, resulting in (J x ) = for Ub = 0. Thus, the origin of the mutual coherence between two molecular 
modes is due solely to two-body coll isions. Furthermore, w e remark that a semiclassical treatment results in (J x ) = 



for all times and all values of Ub /tj l|Mivakawa et al 



2004). Hence, we conclude that the build-up of (J x ) is a purely 



quantum- mechanical effect due to quantum fluctuations. 



The phase distribution of the two wells can be studied using the Pegg-Barnett phase states IjBarnett and Pegg , 
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FIG. 5 Time evolution of P((j>„) for = it, tj = 2.5 and for (a) u b /tj = 0.0032, (b) u b /tj = 0.5623, (c) u b /tj = 56.23. 



199C 


Javanaincn and Ivanov 


1999 


Luis and Sanchez-Soto 


1993 


Pcgg and Barnett 


1988) 



1988). Since the density matrix 



is diagonal in the total number of molecules it is sufficient to consider the relative phase. Figure [S] shows the time 
evolution of the relative phase distribution in three different regimes: (a) Rabi, Ub/tj — 0.0032, (b) Josephson, 
Ub/tj — 0.5623, and (c) Fock Ub/tj = 56.23, for = it, t j = 2.5. Since the vacuum state is taken as the initial state, 
the relative phase at t — is randomly distributed, P{<fi n ) = l/( s + !)■ 

In the Rabi regime, corresponding to Fig|SJa), bimodal phase dis tribution with peaks ar ound both and ±7r builds 



1986). In the Josephson regime, 



up in the characteristic time 7 -1 needed to reach a steady state IjFilipowicz et al 
the relative phase locks around (±7r), for repulsive (attractive) two-body interactions, see Fig. EJV). In contrast to 
these two regimes, in the Fock regime the relative phase distribution becomes almost random for all times, and the 
localized modes in the two wells evolve independently of each other. 

The three regimes of phase distributions correspond to different orders of magnitude of the ratio Ub/tj. The crossover 
of the non-equilibrium steady state from a phase-coherent regime to the random-pha se situation i s reminiscent of the 



super fluid-Mott insulator phase transition for the ground state of an optical lattice l|Fisher et al 



1989: 



Jaksch et al 



1998). Since we consider just two sites, however, there is no sharp transition between these regimes. 



III. PASSAGE TIME STATISTICS OF MOLECULE FORMATION 



We now turn to a second-example that illustrates the understanding of the dynamics of quantum-degenerate atomic 
and molecular systems that can be gained from quantum optics analogies. Here, we consider the first stages of coherent 
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molecular formation via photo-association. Since in such experiments the molecular field is typically in a vacuum 
initially, it is to be intuitively expected that the initial stages of molecule formation will be strongly governed by 
quantum noise, hence subject to large fluctuations. One important way to characterize these fluctuations is in terms 
of the so-called passage time, which is the time it takes to produce, or dissociate, a predetermined number of molecules. 
Quantum noise results in fluctuations in that time, whose probability distribution can therefore be used to probe the 
fluctuations in the formation dynamics. 

Because of the analogy between pairs of fermionic atoms and two-level systems that we already exploited in the 
discussion of the molecular micromaser, one can expect that the problem at hand is somewhat analogous to sponta- 
neous radiation from a sample of two-level atoms, the well-know problem of superradiance. In this section we show 
that this is indeed the case, and use this analogy to study the passage time statistics of molecular formation from 
fermionic atoms. 



A. Model 

We consider again a quantum-degenerate gas of fermionic atoms of mass rrif and spin a —\ , J,, coupled coherently 
to bosonic molecules of mass mj, = 2m t and zero momentum via photo-association. Neglecting collisions between 
fermions and assuming that for short enough times the molecules occupy a single-mode of the bosonic field, this 
system can be described by the boson-fermion model Hamiltonian 



At 



(16) 



where b\ b are molecular bosonic creation and annihilation operators and c ka , Cka are fermionic creation and annihi- 
lation operators describing atoms of momentum hk and spin a. The first and second terms in Eq. (|16fl describe the 
kinetic energy hwk/2 = h 2 k 2 / (2m/) of the atoms and the detuning energy of the molecules respectively, and the third 
term describes the photo-association of p airs of atoms of o pposite momentum into molecules. 



Introducing the pseudo-spin operators ([Anderson . 



1958) analogous to Eq. J7J) 



?(4t £ *T + cl kl c-ki - 1), 



(*fc) f 



c -fc| fcf 



(17) 



the Hamiltonian l|16f) becomes, within an unimportant constant I Barankov and Levitov 



2004: 



Meiser and Mevstrg , 
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2005), 



(18) 



This H amiltonian is known in q uantu m optics as the inhomogeneously broadened (or non-degenerate) Tavis-Cummings 



model l|Tavis and dimming; 



It describes the coupling of an ensemble of two-level atoms to a single- 



mode electromagnetic field. Hence the mapping (|17|) establishes the formal analogy between the problem at 



hand 


and Dicke superradiance, with the caveat that we are dealing with a single bosonic mode ( 


Andreev 


ei^al 


2004 


Barankov and Lev tov . 


2004; 


Javanainen et al. 


2004 


Meiser and Mevstre , 


2005 


Mivakawa and Mevstre 


2005 


Tikhonenkov and Vardi 


200/ 


). Instead of real two-level atoms, pairs of fermionic atoms are now described as effective 



two-level systems whose ground state corresponds to the absence of a pair, \gk) = \0kT,0-ki) and the excited state 
to a pair of atoms of opposite momenta, |efc) = |lfc| , l-fcj,)j in close analogy to the treatment of the atoms in the 
previous section. 

The initial condition consists of the molecular field in the vacuum state and a filled Fermi sea of atoms 



1*7= n #i°>> 
fc<ifc F i 



(19) 



where hp is the Fermi momentum. As such, the problem at hand is in direct analogy to the traditional superradiance 
problem where one starts from an ensemble of excited two-state atoms, as expected from our previous comments. 
Later on we will also consider an initial state containing only molecules and no atoms. This is an important extension 
of the traditional Dicke superradiance system, where the two-level atoms are coupled to all modes of the photon 
vacuum a situation, thereby precluding the possibility of an initial state containing a single, macroscopically occupied 
photon mode unless the system us prepared in a high-Q cavity. 

We assume from now on that the inhomogeneous broadening due to the spread in atomic kinetic energies can be 
ignored. This so-called degenerate approximation is justified provided that the kinetic energies are small compared 
to the atom- molecule coupling energy, (3 = ep/(hx) *C 1, where tp is the Fermi energy. It is the analog of the 
homogeneous broadening limit of quantum optics, and of the Raman-N ath approximation in atomic diffraction. A 



comparison with typical experimental parameters IjHeinzen et 



al 



justified if the number of atoms does not exceed ~ 10 2 — 10 3 IjUvs et al 



2000) shows that the degenerate approximation is 



2005). 



Limiting thus our considerations to small atomic samples, we approximate all Wfc's by up and introduce the collective 
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pseudo-spin operators 



k 



(20) 



obtaining the standard Tavis-Cummings Hamiltonian ( Miv akawa and Mevstre . 



2005; 



Tavis and Cumming; 



1968) 



H = hw F S z + hujbttb + h x (bS + + tfs-). 



(21) 



This Hamiltonian conserves the total spin operator S 2 . The total number of atoms is twice the total spin and hence 
is also a conserved quantity. S z measures the difference in the numbers of atom pairs and molecules. 

Eq. I121H can be diagonalized numerically with reasonable computation times even for relatively large numbers of 
atoms. One can, however, gain significant intuitive insight in the underlying dynamics by finding operator equations 
of motion and then treating the short-time molecular population semiclassically, (hi,) — > rib- To this end we introduce 
the "joint coherence" operators 



f x = (bS+ + tfS-)/2, 

Ty = (bS + -ttS-)/2i, 



and find the Heisenberg equations of motion 



n h 



Tr 



5T„ 



% = ~ST X - x (2S z h b + S+S-^j , 



(22) 

(23) 
(24) 
(25) 



where 6 = ujb — ^>f, so that 2xT x + Shb is a constant of motion. 

In the following, we confine our discussion to the case of S = for simplicity. We thus neglect the contribution of 
T x in Eq. I|25() . In order to better understand the short time dynamics we reexpress S + S~ as 



+ Q- - 



s+s 



-ill + (25 - l)h b + N. 



(26) 



This shows that the operator S + S~ is non- vanishing when the molecular field is in a vacuum and hence can be 
interpreted as a noise operator. Indeed Eqs. (|23|l - (|25|l show that the buildup of the molecular field is triggered only 
by noise if hb = initially. By keeping only the lowest-order terms in hb we can eliminate T y to obtain the differential 
equation 



h b w 2N X 2 (2h b + 1) 



(27) 
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2 3tN 1/2 t 



FIG. 6 Short-time dynamics of (rib}. From left to right, the curves give the l inearized solution 1281 and the full quantum results 



for N = 500, N = 250, and N = 100, respectively. Figure taken from Ref. l|Uvs et al 



2005) 




0.2 ; 



FIG. 7 Passage time distribution for converting 5% of the initial population consisting of only atoms (molecules) into molecules 
(atoms) for N = 500. For initially all atoms: solid line, for initially all molecules: dashed line. 



which, for our initial state, may be solved to yield 



(h b {t)) fa sinli (xVJVt) 



(28) 



Fig. compares the average molecule number (fib) obtained this way, with the full quantum solution obtained by 
direct diagonalization of the Hamiltonian i|21|) for various values of N. The semiclassical approach agrees within 5% 
of the full quantum solution until about 20% of the population of atom pairs has been converted into molecules in all 



cases. 



Next we turn to the passage time statistics. In Fig0we show (solid line) the distribution of times required to produce 
a normalized molecule number n\ ,/N — 0.05 from a sample initially containing N = 500 pairs of atomic fermions, as 
found by direct diagonalization of the Hamiltonian l|21|) . This distribution differs sharply from its counterpart for the 
reverse process of photodissociation from a molecular condensate into fermionic atom pairs, which is plotted as the 
dashed line in Fig. [7| In contrast to photo-association, this latter process suffers significantly reduced fluctuations. 

To understand the physical mechanism leading to this reduction in fluctuations we again turn to our short time 



17 




FIG. 8 Effective potential for a system with TV ^> 1. The circle (square) corresponds to a n initial state wit h all fermionic atoms 



(molecules). The part of the potential for rib < is unphysical. Figure taken from Ref. IjUvs et al 



2005) 



semi-classical model. Within t his approximation, the Heisenb erg equations of motion l|23(l - l|24(l can be recast in the 



form of a Newtonian equation IjMivakawa and Mevstrg , 



2005) 



d 2 n b 



dU{n b ) 



(29) 



dt 2 drib ' 

where the cubic effective potential U(rib) is plotted in Fig. |H1 (Note we have now kept all orders in n^.) In case 
the system is initially composed solely of fermionic atoms, n&(0) = 0, the initial state is dynamically unstable, with 
fluctuations having a large impact on the build-up of rij,. In contrast, when it consists initially solely of molecules, 
715 = N, the initial state is far from the point of unstable equilibrium, and n b simply "rolls down" the potential in a 
manner largely insensitive to quantum fluctuations. This is a consequence of the fact that the bosonic initial state 
provides a mean field that is more amenable to a classical description. Hence, while the early stages of molecular dimer 
formation from fermionic atoms are characterized by large fluctuations in formation times that reflect the quantum 
fluctuations in the initial atomic state, the reverse process of dissociation of a condensate of molecular dimers is largely 
deterministic. The diminished fluctuations in this reversed process is peculiar to the atom-molecule system and not 
normally considered in the quantum optics analog of Dicke superradiance. 



IV. COUNTING STATISTICS OF MOLECULAR FIELDS 

An important quantum mechanical characteristic of a quantum field is its counting (or number) statistics. In this 
section we show how the similarity of the coherent molecule formation with quantum optical sum-frequency generation 
can be used to determine the counting statistics of the molecular field. In particular we show how the counting statistics 
depends on the statistics of the atoms from which the molecules are formed. Besides being interesting in its own right, 
such an analysis is crucial for an understanding of several recent experiments that used a "projection" onto molecules 



to detect BCS superfluidity in fermionic systems l)Regal et al 



2004a; 



Zwierlein et al 
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2004). Our work shows that 



the statistical properties of the resulting molecular field indeed reflect properties of the initial atomic state and are a 
sensitive probe for superfluidity. 

As before, we restrict our discussion to a simple model in which all the molecules are generated in a single mode. We 
use time dependent perturbation theory to calculate the number of molecules formed after some time t, n{t), as well 
as the equal-time second-order correlation function g( 2 >(t,t). We also integrate the Schrodinger equation numerically 
for small numbers of atoms, which allows us to calculate the complete counting statistics P n . 



A. BEC 



Consider first a cloud of weakly interacting bosons well below the condensation temperature T c . It is a good 
approximation to assume that all atoms are in the condensate, described by the condensate wave function fx ) . The 



coupled system of atoms and molecules is described by the effective two-mode Hamiltonian IjAnglin and Vardi , 



Javanainen and Mackie , 



1999) 



2001 



Hbec = hSfcb + h X (b^c 2 + Sc t2 ) 



(30) 



where b, w and c, & are the bosonic annihilation and creation operators for the molecules and for the atoms in the 
condensate, respectively, 8 is the detuning between the molecular and atomic level, and H\ is the effective coupling 
constant. 

Typical experiments start out with all atoms in the condensate and no molecules, corresponding to the initial state, 



m=o)) = ^=\o). (si) 

where N a = 2A< r max is the number of atoms, N max is the maximum possible number of molecules and |0) is the vacuum 
of both molecules and atoms. We can numerically solve the Schrodinger equation for this problem in a number basis 
and from that solution we can determine the molecule statistics P n (t) . The results of such a simulation are illustrated 
in Fig. |SJ which shows P n (t) for 30 initial atom pairs and 6 = 0. Starting in the state with zero molecules, a 
wave-packet-like structure forms and propagates in the direction of increasing n. Near iV max the molecules begin to 
dissociate back into atom pairs. 

We can gain some analytical insight into the short-time dynamics of molecule formation by using first-order per- 



turbation theory l)Kozierowski and R. Tanas , 



1977 



Mande 



1982). We find for the mean molecule number 



n(t) = (xt) 2 2N max (2N max - 1) + 0{{ X t) 2 ) 



and for the second factorial moment 



+ 0(N t 



2 ) 

max / 
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(32) 



(33) 



For iV max large enough we have g^ 2 '(ti,t2) — ► 1, the value characteristic of a Glauber coherent field. From g( 2 > and 
n{t) we also find the relative width of the molecule number distribution as 



1. 



(34) 



It approaches n 1 I' 2 in the limit of large -/V max , typical of a Poisson distribution. This confirms that for short enough 
times, the molecular field is coherent in the sense of quantum optics. 



B. Normal Fermi gas 

We now turn to the case of photo-association from two different species of non-interacting ultra-cold fermions. The 
two species are again denoted by spin up and down. At T = 0, the atoms fill a Fermi sea up to an energy /i. Weak 
repulsive interactions give rise only to minor quantitati ve modifications tha t we ignore. We refer to this system of 



non-interacting Fermions as a normal Fermi gas (NFG) IjLandau et al 



1980). 



As before we assume that atom pairs are coupled only to a single mode of the molecular field, which we assume 
to have zero momentum for simplicity. Then, using the mapping to pseudo spins Eq. (|17fl we find that the system 
is again described by the inhomogeneously broadened Tavis-Cummings Hamiltonian Eq. I|18fl . However, in contrast 
to the previous case, we do not assume that the fermionic energies are approximately degenerate, in order to be able 
compare the results to the BCS case, where the kinetic energies are essential. 

Figure^]shows the molecule statistics obtained this way. The result is clearly both qualitatively and quantitatively 
very different from the case of molecule formation from an atomic BEC. From the Tavis-Cummings model analogy 
we expect that for short times the statistics of the molecular field should be chaotic, or "thermal", much like those 
of a single-mode chaotic light field. This is because each individual atom pair "emits" a molecule independently and 
without any phase relation with other pairs. That this is the case is illustrated in the inset of Fig. I1UI which fits the 
molecule statistics at selected short times with chaotic distributions of the form 



Pn 



D -n/ (n) 



thermal 



V e -n/{n) 



(35) 
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FIG. 9 Number statistics of molecules formed from a BEC with iV max = 30 and S = 0. 

The increasing 'pseudo-temperature' (n) corresponds to the growing average number of molecules as a function of 
time. 

As before we determine the short-time properties of the molecular field in first-order perturbation theory. We find 
for the mean number of molecules 

n(t) = (xt) 2 2N a . (36) 

It is proportional to N a , in contrast to the BEC result, where n was proportional to N%, see Eq. (J22J- This is 
another manifestation of the independence of all the atom pairs from each other: While in the BEC case the molecule 
production is a collective effect with contributions from all possible atom pairs adding constructively, there is no such 
collective enhancement in the case of Fermions. Each atom can pair up with only one other atom to form a molecule. 
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FIG. 10 Number statistics of molecules formed from a normal Fermi gas. This simulation is for N a = 20 atoms, the detuning 
is 8 = 0, the Fermi energy is fi = O.lhx and the momentum of the i-th pair is \ki\ = (i — 1)2/cf /(N a /2 — 1). The inset shows 
fits of the number statistics to thermal distributions for various times as marked by the thick lines in the main figure. 

For the second factorial moment wc find 



5 ( 2 )(^,< 2 ) = 2(1 



2N a 



(37) 



which is close to two, typical of a chaotic or thermal field. 



C. Fermi gas with superfluid component 

Unlike repulsive interactions, attractive interactions between fermions have a profound impact on molecule forma- 
tion. It is known that such interactions give rise to a Cooper instability that leads to pairing and drastically changes 
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FIG. 11 Number statistics of molecules formed from a Fermi gas with pairing correlations. For this simulation the detuning is 
5 — 0, the Fermi energy is /i = O.lg and the background scattering strength is V = 0.03x resulting in N a ~ 9.4 atoms and a 
gap of A ~ 0.15x- The momenta of the atom pairs are distributed as before in the normal Fermi gas case. 

the qualitative properties of the atomic system. The BCS reduced Hamiltonian is essentially the inhomogeneously 
broadened Ta vis-Cumming s Hamiltonian Ijl8|l with an additional term accounting for the attractive interactions be- 



tween atoms (Kittc 



1987) 



H = r ^kK + M r b + h x Y, (ttt k + ba+) fit** ■ 



(38) 



fc.fc' 



The approximate mean-field ground state \BCS) is found by minimizing (Hbcs ~ (*N) in the standard way. The 
dynamics is then obtained by numerically integrating the Schrodinger equation with |BCS) as the initial atomic state 
and the molecular field in the vacuum state. 

Figure ITT1 shows the resulting molecule statistics for V — O.QShx, which corresponds to a gap of A = O.lbhx for the 
system at hand. Clearly, the molecule production is much more efficient than it was in the case of a normal Fermi 
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gas. The molecules are produced at a higher rate and the maximum number of molecules is larger. The evolution of 
the number statistics is reminiscent of the BEC case. This also shows that the qualitative differences seen between 
the normal Fermi gas and a BEC in the previous section cannot be attributed to inhomogeneous broadening and the 
resulting dephasing alone but are instead a result of the different coherence properties of the atoms. 
The short-time dynamics is again obtained in first-order perturbation theory, which gives now 



In addition to the term proportional to N a representing the incoherent contribution from the individual atom pairs 
that was already present in the normal Fermi gas, there is now an additional contribution proportional to (A/V) 2 . 
Since (A/V) can be interpreted as the number of Cooper pairs in the quantum-degenerate Fermi gas, this term 
can be understood as resulting from the coherent conversion of Cooper pairs into molecules in a collective fashion 
similar to the BEC case. The coherent contribution results naturally from the nonlinear coupling of the atomic field 
to the molecular field. This nonlinear coupling links higher-order correlations of the molecular field to lower-order 
correlations of the atomic field. For the parameters of Fig. EI A/V « 6.5 so that the coherent contribution from the 
Cooper pairs clearly dominates over the incoherent contribution from the unpaired fermions. Note that no signature 
of that term can be found in the momentum distribution of the atoms themselves. Their momentum distribution is 
very similar to that of a normal Fermi gas. The short-time value of ff^(*i>*2)> shown in Fig. ^1 decreases from 
the value of Eq. (|37J) for a normal Fermi gas at A = down to one as A increases, underlining the transition from 
incoherent to coherent molecule production. 

V. MOLECULES AS PROBES OF SPATIAL CORRELATIONS 

The single-mode description of the molecular field of the previous section results in the loss of all information about 
the spatial structure of the atomic state. In this final section we adopt a complementary view and study the coupled 
atom-molecule system including all modes of the molecular and atomic field so as to resolve their spatial structure. 
This problem is too complex to admit an exact solution, hence we rely entirely on perturbation theory. 

One of the motivations for such studies is the on-going experimental efforts to study the so-called BEC-BCS 
crossover. A difficulty of these studies has been that they necessitate the measurement of higher-order correlations of 
the atomic system. While the momentum distribution of a gas of bosons provides a clear signature of the presence of a 
Bose-Einstein condensate, the Cooper pairing between fermionic atoms in a BCS state hardly changes the momentum 




(39) 
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FIG. 12 g (2) (0+,0 + ) as a function of the gap parameter A. (Figure taken from Ref. IjMeiser and MevstreL 12005")) 



distribution or spatial profile as compared to a normal Fermi gas. This poses a significant experimental challenge, 
since the primary techniques for probing the state of an ultra-cold gas are either optical absorption or phase contrast 
imaging, which directly measure the spatial density or momentum distribution following ballistic expansion of the 
gas. In the strongly interacting regime very close to the Feshbach resonance, evidence for fermionic sup erfluidity was 



obtained by projecting the a tom pairs onto a molecular state by a rapid sweep through the resonance l|Regal et al 



2004b; 



Zwierlein et al 



rf sp e ctroscopy (|C 



2004; 



Kinast et al 



n et al 



2004) . Mor e direct evidence of the gap in the excitation spectra due to pairi ng was obtained by 



2004) and by measurements of the collective excitation frequencies l|Bartenstein et al 



2004). Finally, the superfluidity of ultra-cold fermions i n the strongly interact ing regime has 



recently been impressively demonstrated via the generation of atomic vortices l|Zwierlein et al 



2005). 
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Still, the detection of fermionic superfluidity in the weakly interacting BCS regime remains a challenge. The di- 
rect detection of Cooper pairing requires the measurement of second-order or higher atomic correlation functions. 
Several researchers have proposed and implemen t ed schemes that a llow one to measure higher order corre l ations 



Altmanet 



Rega l et al 



al 



2004; 
13) 1 



Bach and Rzazcwski 



2004; 



Burt et al 



1997 



Cacciapuoti et al 



2003; 



Hellweg et al 



2003; 



2004b) but those methods are still very difficult to realize experimentally. While the measurement 



of higher-order correlations is challengi ng already for bosons, the theory of these correlation s has been established a 



despite some efforts ([Cahill and Glauber 



long time ago by Gla uber for photons IIGlauber 



1963a b; 



Naraschewski and Glauber 



1999). For fermions however, 



1999) a satisfactory coherence theory is still missing. 



From the previous section we know that one can circumvent these difficulties by making use of the nonlinear coupling 
of atoms to a molecular field. The nonlincarity of the coupling links first-order correlations of the molecules to second- 
order correlations of the atoms. Furthermore the molecules are always bosonic so that the well-known coherence theory 
for bosonic fields can be used to characterize them. Considering a simplified model with only one molecular mode, it 
was found that the molecules created that way can indeed be used as a diagnostic tool for second-order correlations 
of the original atomic field. 

We consider the limiting case of strong atom-molecule coupling as compared to the relevant atomic energies. The 
molecule formation from a Bose-Einstein condensate (BEC) serves as a reference system. There we can rather easily 
study the contributions to the molecular signal from the condensed fraction as well as from thermal and quantum 
fluctuations above the condensate. The cases of a normal Fermi gas and a BCS superfluid Fermi system are then 
compared with it. We show that the molecule formation from a normal Fermi gas and from the unpaired fraction of 
atoms in a BCS state has very similar properties to those of the molecule formation from the non-condensed atoms 
in the BEC case. The state of the molecular field formed from the pairing field in the BCS state on the other hand 
is similar to that resulting from the condensed fraction in the BEC case. The qualitative information gained by the 
analogies with the BEC case help us gain a physical understanding of the molecule formation in the BCS case where 
direct calculations are difficult and not nearly as transparent. 



A. Model 

We consider again the three cases where the atoms are bosonic and initially form a BEC, or consist of two species 
of ultra-cold fermions (labeled again by a =|, I), with or without superfluid component. In the following we describe 
explicitly the situation for fermions, the bosonic case being obtained from it by omitting the spin indices and by 
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replacing the Fermi field operators by bosonic field operators. 

Since we are primarily interested in how much can be learned about the second-order correlations of the initial 
atomic cloud from the final molecular state, we keep the physics of the atoms themselves as well as the coupling to 



1 Chiofalo et al 



2002 


Holland ct al. 


2001 


Timmermans et al. 


1999) 



& = ET £ ^ 5k - + E^4«k + ^ 1/2 E ^( k 2 - k 04 2 <A lc 



k.tr k ki ,k-2 ,<r 

+ 2V E £ k 1+qT 5 k 2 -qiCk 2 ic klT + % I E aqC q /2+kiC q/2 _ kT + H.c. (40) 

q,ki,k 2 \q,k ) 

The kinetic energies Uk are defined as before, V is the quantization volume, Z7t r (k) = V^ 1 ^ 2 j v G? 3 ire _lkx [/tr(x) is the 
Fourier transform of the trapping potential £/tr(r) and Uq = Aith 2 a/mf is the background scattering strength with o 
the background scattering length. The coupling constant g between atoms and molecules is, up to dimensions, equal 
to x °f the previous sections. 

We assume that the trapping potential and background scattering are relevant only for the preparation of the 
initial state before the coupling to the molecules is switched on at t = and can be neglected in the calculation of 
the subsequent dynamics. This is justified if Hg^/N 3> Uon, huji where n is the atomic density, N the number of 
atoms, and uji are the oscillator frequencies of the atoms in the potential J7 tr (r) that is assumed to be harmonic. 
Experimentally, the interaction between the atoms can effectively be switched off by ramping the magnetic field to a 
position where the scattering length is zero, so that this assumption is fulfilled. 

Regarding the strength of the coupling constant g, two cases are possible: hgy/N can be much larger or much smaller 
than the characteristic kinetic energies involved. For fermions the terms broad and narrow resonance have been coined 
for the two cases, respectively, and we will use these for bosons as well. Both situations can be realized experimentally, 
and they give rise to different effects. For strong coupling the conversion process needs not satisfy energy conservation 
because of the energy time uncertainty relation. For weak coupling energy conservation is enforced. This energy 
selectivity can be useful in certain situations because it allows one to resolve additional structures in the atomic state. 
The analysis of this case is fairly technical, ho wever. Therefore we only consider the case of strong coupling and refer 



the interested reader to (Meiser et al 



2005bJ) for details of the calculations and the case of weak coupling. 
First-order time-dependent perturbation theory requires that the state of the atoms does not change significantly 
and consequently, only a small fraction of the atoms are converted into molecules. It is reasonable to assume that 



27 

this is true for short interaction times or weak enough coupling. Apart from making the system tractable by analytic 
methods there is also a deeper reason why the coupling should be weak: Since we ultimately wish to get information 
about the atomic state, it should not be modified too much by the measurement itself, i.e. the coupling to the 
molecular field. Our treatment therefore follows the same spirit as Glauber's original theory of photon detection, 
where it is assumed that the light-matter coupling is weak enough that the detector photocurrent can be calculated 
using Fermi's Golden rule. 



B. BEC 



We consider first the case where the initial atomic state is a BEC in a spherically symmetric harmonic trap. We 
assume that the temperature is well below the BEC transition temperature and that the interactions between the 
atoms are not too strong. Then the atomic system is described by the field operator 



Vj(x) = V'o(x)c + (5^(x), 



(41) 



where ipo ( x ) is the condensate wave function and c is the annihilation operator for an atom in the condensate. In 
accordance with the assumption of low temperatures and weak interactions we do not distinguish between the total 
number of atoms and the number of atoms in the condensate. The fluctuations 5ip(x) are small and those with 



comp arable to Rtf can be neglected l|Bergeman et al 



1999) 



2000; 


D. A. W. Hutchinson and Zarcmba 


1997 


Rcidl et al. 



We are interested in the momentum distribution of the molecules 

n(p,t) = (bl(t)b p (t)) 



(42) 



which for short times, t, can be calculated using perturbation theory. In the broad resonance limit we ignore the 
kinetic energies and find 



n BEC (p, t) = (gt) 2 N(N - l)V iP 2 (p) + (gtfAN \ — (bc% (x)5c p (x)) 



dh 



(43) 



where the expectation value in the last term includes a thermal average. From this expression we see that our approach 
is justified if (\/~N gt) 2 <C 1 because for such times the initial atomic state can be assumed to remain undepleted. The 
first term in Eq. (|43|) is the contribution from condensed atoms and the second term comes from uncondensed atoms 



2s 



above the condensate. The contribution from the condensate can be evaluated in closed form in the Thomas Fermi 
approximation for a spherical trap. The contribution from the therm al atoms can be calcu lated using the local density 



approximation. The details of this calculation can be found in Ref. l|Meiser et al 



2005b) 



The momentum distribution (|43[) is illustrated in Fig. The contribution from the condensate is a collective effect, 
as indicated by its quadratic scaling with the atom number. It clearly dominates over the incoherent contribution 
from the fluctuations, which is proportional to the number of atoms and only visible in the inset. The momentum 
width of the contribution from the condensate is roughly /i27r/i?TF which is much narrower than the contribution 
from the fluctuations, whose momentum distribution has a typical width of where £ = (87ran) _1 ^ 2 is the healing 
length. This is a case where coherence properties of the atoms can be read off the momentum distribution of the 
molecules: The narrow momentum distribution of the molecules is only possible if the atoms were coherent over 
distances ~ Rtf- At this point this is a fairly trivial observation and the same information could have been gained 
by looking directly at the momentum distributi on of the atoms, which i s after all how Bose-Einstein condensati on was 



detected already in the very first experiments (And erson et al 



1995; 



Bradley et al 



1995; 



Davis et al 



Still 



we mention it because it will be very interesting (indeed interesting enough to motivate this whole work!) to contrast 
this situation to the BCS case below. 

Using the same approximation scheme we can calculate the second-order correlation. If we neglect fluctuations we 
find 



ffBEcCPi'^Ps^) = 1 



N 



0{N~ 2 ). 



(44) 



For N — > oo this is very close to 1, which is characteristic of a coherent state. This result implies that the number 
fluctuations of the molecules are very nearly Poissonian. The fluctuations lead to a larger value of , making the 
molecular field partially coherent, but their effect is only of order 0(N^ 1 ). 



C. Normal Fermi gas 



We treat the gas in the local density approximation where the atoms locally fill a Fermi sea 



\NFG) = J] 4|0) 

|k|<fcjr(x) 



(45) 



with local Fermi momentum /tfcp(x) and |0) being the atomic vacuum . It is related to the local density of the atoms 



in the usual way l|D. A. Butts and D. S. Rokshar 



1997 



Landau et al 



1980) 
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FIG. 13 Momentum distribution of molecules formed from a BEC (red dashed line) with a — 0.1a osc and T — 0.1T C and a BCS 
type state with fc^a = 0.5 and a osc = 5fc^ 1 (0) (blue solid line), both for N = 10 5 atoms. The BCS curve has been scaled up by 
a factor of 20 for easier comparison. The inset shows the noise contribution for BEC (red dashed) and BCS (blue) case. The 
latter is simply the momentum distribution of molecules formed from a normal Fermi gas. The local density approximation 
treatment of the noise contribution in the BEC case is not valid for momenta smaller than 2-7r/£ (indicated by the red dotted 
line in the inset). Note that the coherent contribution is larger than the noise contribution by five orders of magnitude in the 
BEC case and three orders of magnitude in the BCS case. 

The momentum distribution and second order correlation function are readily found in perturbation theory. The 
momentum distribution is shown in the inset in Fig. 1131 The total number of molecules scales only linear with the 
number of atoms, meaning that, in contrast to the BEC case, the molecule formation is non-collective. Each atom 
pair is converted into a molecule independently of all the others and there is no collective enhancement. Furthermore 

1/3 

the momentum distribution of the atoms is much wider than in the BEC case. It's width is of the order of hn 
indicating that the atoms are correlated only over distances comparable to the inter atomic distance. 
Similarly, we find for the local value of at position x 

Cc(P> x >*) = Sioc(P>*;P'*> x ) 

-'('-JsbO- (46) 

where N e g is the number of atoms that are allowed to form a molecule on the basis of momentum conservation. For 
large N e s approaches 2 which is characteristic of a thermal field. Indeed, using the analogy with an ensemble of 
two level atoms coupled to every mode of the molecular field provided by the Tavis-Cummings model, it is easy to 
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see that the entire counting statistics is thermal. 



D. BCS state 

Let us finally consider a system of fermions with attractive interactions, Uq < 0, at temperatures well below the 
BCS critical temperature. As is well known, for these temperatures the attractive interactions give rise to correlations 
between pairs of atoms in time reversed states known as Cooper pairs. We assume that the spherically symmetric 
trapping potential is sufficiently slowly varying that the gas can be treated in the local density approximation. More 
quantitatively, the local density approximation is valid if the size of the Cooper pairs, given by the correlation length 

A(r) = u_F(r)/7rA(r), 

is much smaller than the oscillator length for the trap. Here, vp(r) is the velocity of atoms at the Fermi surface and 
A(r) is the pairing field at distance r from the origin, which we take at the center of the trap. Loosely speaking, in the 
local density approximation the ground state of the atoms is determined by repeating the variational BCS-calculation 
of the p revious section in a sm all volume at every position x. A thorough discussion of this calculation can be found 



in Rcf. (Houbiers et al 



1997). 

We find the momentum distribution of the molecules from the BCS-type state by repeating the calculation done 
in the case of a normal Fermi gas. For the BCS wave function, the relevant atomic expectation values factorize into 
normal and anomalous correlations. The normal terms are proportional to densities and are already present in the 
case of a normal Fermi gas while the anomalous contributions are proportional to the gap parameter. The momentum 



distribution of the molecules becomes 



^2 ( fi p/2+k,lCp/2-k,T 



•n NFG (p,i). (47) 



k 

The first term is easily shown to be proportional to the square of the Fourier transform of the gap parameter. 
Since the gap parameter is slowly varying over the size of the atomic cloud, this contribution has a width of the order 
of H/Rtf, in complete analogy with the BEC case above. The total number of atoms in the first contribution is 
proportional to the square of the number of Cooper pairs, which is a macroscopic fraction of the total atom number 
well below the BCS transition temperature. That means that this contribution is a collective effect. The second term 
is the wide and incoherent non-collective contribution already present in the case of a normal Fermi gas. It is very 
similar to the thermal noise in the BEC case as far as its coherence properties are concerned. 
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For weak interactions such that the coherent contribution is small compared to the incoherent contribution, the 
second order correlations are close to those of a normal Fermi gas given by Eq. (|46Jl . g^ 2 '(p, x, t) k, 2. However, in the 
strongly interacting regime, fci?H ~ 1, and large N , the coherent contribution from the paired atoms dominates over 
the incoherent contribution from unpaired atoms. In this limit one finds that the second-order correlation is close 
to that of the BEC, g( 2 '(p,x, t) w 1. The physical reason for this is that at the level of even order correlations the 
pairing field behaves just like the mean field of the condensate. This is clear from the factorization property of the 
atomic correlation functions in terms of the normal component of the density and the anomalous density contribution 
due to the mean field. In this case, the leading order terms in N are given by the anomalous averages. 

To summarize, molecules produced from an atomic BEC show a rather narrow momentum distribution that is 
comparable to the zero-point momentum width of the BEC from which they are formed. The molecule production 
is a collective effect with contributions from all atom pairs adding up constructively, as indicated by the quadratic 
scaling of the number of molecules with the number of atoms. Each mode of the resulting molecular field is to a very 
good approximation coherent (up to terms of order 0(1/N)). The effects of noise, both due to finite temperatures 
and to vacuum fluctuations, are of relative order 0(1/N). They slightly increase the g^ and cause the molecular 
field in each momentum state to be only partially coherent. 

In contrast, the momentum distribution of molecules formed from a normal Fermi gas is much broader with a 
typical width given by the Fermi momentum of the initial atomic cloud, i.e. the atoms are only correlated over an 
interatomic distance. The molecule production is not collective as the number of molecules only scales like the number 
of atoms rather than the square. In this case, the second-order correlations of the molecules exhibit super-Poissonian 
fluctuations, and the molecules are well characterized by a thermal field. 

The case where molecules are produced from paired atoms in a BCS-like state shares many properties with the 
BEC case: The molecule formation rate is collective, their momentum distribution is very narrow, corresponding to a 
coherence length of order Rtf, and the molecular field is essentially coherent. The non-collective contribution from 
unpaired atoms has a momentum distribution very similar to that of the thermal fluctuations in the BEC case. 

VI. CONCLUSION 

In this paper we have used three examples to illustrate the profound impact of quantum optics paradigms, tools 
and techniques, on the study of low-density, quantum-degenerate atomic and molecular systems. There is little doubt 
that the remarkably fast progress witnessed by that field results in no little part from the experimental and theoretical 
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methods developed in quantum optics over the last decades. It is therefore fitting, on the occasion of Herbert Walther's 
seventieth birthday, to reflect on the profound impact of the field that he has helped invent, and where he has been 
and remains so influential, on some of the most exciting developments in AMO science. 
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